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Abstract 



It is shown that the Abrikosov's vortex solution or its corresponding two- 
particle pair potential is not the solution of the self-consistency equation in 
Gor'kov's formalism. Since the self-consistency equation leads to a superpo- 
sition of different types of off-diagonal long-range order (ODLRO) instead of 
one type of ODLRO only, it may not handle the vortex problem appropriately. 
A possible resolution is suggested. 



PACS numbers: 74.50.+r, 74.60.-w, 74.60.Ec 
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I. INTRODUCTION 



Recent STM experiments 1 " 3 show that the microscopic structures of the vortex state 
in type II and high T c superconductors are much more complicated than expected. In 
particular, the vortex core and the lattice structure depend strongly on the microscopic 
details of the samples. For instance, using NbSe 2 samples Hess and his collaborators 1,2 
found unexpected zero-bias peaks and sixfold star-shaped structures in the vortex core states. 
These experiments revived much theoretical interest. 4 " 6 However, a satisfactory quantitative 
explanation is still not available. 7 Maggio-Aprile et al. 3 also reported that the vortices in 
Y BdiCuzOis arrange in an oblique lattice, which remains unexplained. On theoretical 
side, Thouless et al. 8 noted that there is still a lot of controversy in the dynamics of vortices. 

It seems that there is some fundamental problem in our theoretical understanding of the 
vortex state in superconductors. In this paper, we point out that Gor'kov's Green's function 
formalism or the Bogoliubov-de Gennes equations may not handle the vortex problem prop- 
erly, since the self-consistency equation gives rise to an incoherent superposition of different 
types of off-diagonal long-range-order (ODLRO) instead of a coherent vortex state. The 
same problem was already found in Gor'kov and Galitskii's (GG) 9 solution for the d-wave 
pairing BCS theory. GG allowed a superposition of several distinct types of the off-diagonal 
long-range-order, which was proven to be invalid. 10 " 14 As Anderson noted, 13 this problem 
occurs in many cases even in discussing flux quantization. In fact, it also hindered correct 
understanding of the effects of magnetic impurities 15,16 and weak localization 17 " 20 on super- 
conductivity. Balian, Nosanow, and Werthamer 11 showed that this problem is caused by 
the difficulty in the method of Green's functions for a manybody system. It is inspiring to 
remind their statement: Thus the Green's function method as usually formulated is not a 
complete dynamical description of the system, and requires in addition, some criterion to 
distinguish these extraneous solutions from the correct one. The criterion may be provided 
by the physical constraint of the Anomalous Green's function, which is nothing but the 
pairing constraint. 17,18 
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The vortex problem may be very complicated; so we use the following simplifications 
which are not crucial in our discussion though. 

(i) We consider only the lowest Landau Level in the presence of the very high magnetic 
fields, i.e., 



where uj c and wo are the cyclotron frequency and the Debye frequency, respectively. 

(ii) The Zeeman splitting is disregarded. 

(iii) We use the Landau gauge where 



The z-axis motion will be suppressed. 

(iv) We consider the self-consistency equation near T c . 

II. LANDAU LEVELS AND ABRIKOSOV'S VORTEX SOLUTION 

We consider a rectangular sample with sides L x and L y . For the Hamiltonian 

2 / _ eBx\2 
}J — JjL _|_ c I 

2m 2m 

the eigenfunctions are given by 



LO c > UJ D , 



(1) 



A x — ; Ay — Bx. 



(2) 




(4) 



where 




(5) 



and 



Here N Q is a normalization constant and we have used the periodic boundary condition along 
the y-direction. If the x dimensions of the system are confined to —L x /2 < x < L x /2, n is 
determined by the condition 

" Y < nqf < ^. (7) 

(We neglect the effect of boundary on N Q .) 
The famous Abrikosov's vortex solution is 21 

oo 

*(x,y) = £ C n e^yex P [-(x-qn£ 2 ) 2 /2£% (8) 

n=— oo 

with the periodicity conditions, C n+ i = C n for a square lattice and C n+2 = C n for a tri- 
angular lattice, 22 respectively. The C n and t are adjusted to minimize the free energy in 
the Ginzburg-Landau theory. Note that the solution is a linear combination of one-particle 
eigenfunctions. 

III. SELF-CONSISTENCY EQUATION 

In terms of the Lowest Landau Level wavefunctions, the normal state Green's function 
G is written as 

G w (l,2) = V ^j^(^ > (9) 

where e Q = \Tiuj c — /i. 23 /x is the Fermi energy and uj are the Matsubara frequencies. The 
Gor'kov's self-consistency equation is then given by 24 

A( Xl , yi ) — VT [ Y, Gl(l,2)Gi w (l, 2) A(x 2 ,y 2 )dx 2 dy 2 . (10) 

According to the Gor'kov's microscopic derivation of the Ginzburg-Landau equations, 24 
the order parameter \l/(r) is proportional to the (two-body) pair potential A(r): 



*(r) = A(rW7C(3)A47rT c , (11) 
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where ( is Riemann's zeta function and N is the electron number density in the normal 
metal. Consequently, the Abrikosov's vortex solution Eq. (8) should be the solution of the 
above self-consistency equation. Unfortunately, this is not the case. For example, if we 
substitute the pair potential 

A(x, y) ~ C ie iqy exp[-{x - q£ 2 ) 2 /2£ 2 }, (12) 

into the self-consistency equation, we obtain the different form of the pair potential 

A(x,y) ~ C ie tqy exp[-(x - q£ 2 /2) 2 /£ 2 ]. (13) 

In other words, the self-consistency condition is not satisfied. This evidence cast serious 
doubt on the Gor'kov's microscopic derivation of the Ginzburg-Landau equations. In fact, 
this difficulty is anticipated since the Abrikosov's vortex solution is a linear combination 
of one-particle eigenfunctions while the pair potential consists of the multiplication of two 
one-particle eigenfunctions. 

IV. ABSENCE OF THE VORTEX SOLUTION 

Another possible form of the pair potential is 

A(x,y) = J2&ne iqny exp[-(x - qn£ 2 /2) 2 jt 2 \. (14) 

n 

Now we show that this pair potential is not the solution of the self-consistency equation, 
either. First, we consider the constant pair potential in the y-direction, i.e., 

A(x, y) ~ A exp(-x 2 /£ 2 ). (15) 

Upon substitution of this simple solution into Eq. (10), we find pairing between n j and 
0_„ I and the transition temperature T c is determined by 

at2 -2n 2 q 2 £ 2 
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We could have obtained the same equation from the BCS theory with the pairing matrix 
elements, 

Vnw = VJ 0;(r)C n (r)0_ n ,(r)^(r)dr 
The manybody ground state is, then 

0BCS = nK + M^T,0-nl)]|O> (18) 

n 

Second, we consider the pair potential for n — 1 in Eq. (14), 

A(x,y) ~ A ie tqy exp[-(x - q£ 2 /2) 2 /£ 2 ). (19) 

It is straightforward to show that this pair potential leads to pairing between <p n j and 
0-n+i I an d the resulting transition temperature T' c is determined by 

N 2 -2(n-l/2)V^ 2 

— ,20) 

In the BCS theory, the corresponding pairing matrix elements are 

Km' = V J 0;(r)01 n+1 (r)0„ n / +1 (r)0 n /(r)rfr 

= yi^L e -(n-l/2) 2 q 2 ^ 2 e (n'-l/2) 2 g^ 2 ^1) 

V2 ' 

Note that the transition temperatures may be extremely small due to the exponential factors 
and T c and T' c are different. The ground state is now 

&BCS = IlK + ^(0n T,0-n+l l)]|0 > . (22) 

n 

If we combine the two solutions, we obtain the pair potential 

A(x,y) ~ A exp(-x 2 /£ 2 ) + A ie tqy exp[-(x - qf/2) 2 /£ 2 }. (23) 
Inserting Eq. (23) into Eq. (10) one finds 
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Aoexpi-xl/f) =VT f J2Gl(l,2)GU^^)^xp(-xl/£ 2 )dx 2 dy 2 , (24) 

and 

A 1 e^exp[-(x 1 -q£ 2 /2) 2 /£ 2 ] = VT f ]T G£(l, 2)G^(1, 2) 

x A ie iW2 ea;p[-(a; 2 - q£ 2 /2) 2 /£ 2 ]dx 2 dy 2 , (25) 

since the two solutions are linearly independent. Notice that we can not find the tempera- 
ture at which the two different types of the condensation occur simultaneously due to the 
difference in the pairing matrix elements as shown above. Nevertheless, as Galitskii 25 sug- 
gested, it is tempting to write the resulting manybody state as a combination of the above 
ground states, that is, 

&BCS + 4>bcs = 

UK + V n {<t> n T, 4>-n l)]|0 > + IIK + M<f>n T, <f>-n+l I)] |o > . (26) 

n n 

This combination is just the so-called incoherent superposition of different types of off- 
diagonal long-range-order(ODLRO), 9-14 which does not correspond to any real physical 
state. Hone 10 actually demonstrated the impossibility of constructing a complete hierarchy 
of Green's functions in such a case. 
For the pair potential 

A(x,y) ~ A 2 e t2qy exp[-(x - q£ 2 ) 2 /£ 2 }, (27) 

one finds the T c equation 



/V 2 _ ~2(n-l) 2 qH 2 

1=VT -viS?^r (28) 

When the electrons are confined in the x-direction as assumed here, the T c is different from 
those of the previous cases. Whereas for an infinite system or the system with the periodic 
boundary condition in the x-direction, the T c may be the same as that for the pair potential 
corresponding to n = in Eq. (14). We have then T c for A , A 2 , A 4 , • • • and T c '(^ T c ) for 
Ai, A 3 , A 5 , • ■ •. It is interesting to note that the Abrikosov solution with C n+2 = C n has a 
lower energy than that with C n+ i = C n . 



Thus, the self-consistency equation does not allow the Abrikosov's vortex solution or 
the corresponding two-particle pair potential Eq. (14). It is obvious that adding different 
types of off-diagonal long-range-order (ODLRO) in Gor'kov's formalism does not lead to a 
coherent vortex state. On the other hand, previous workers first sum over the eigenstate (f> n 
in the one-particle Green's function and then consider the self-consistency condition, which 
fails to take into account the two-particle correlations correctly. 



V. DISCUSSION 

To describe the vortex state, we need to devise a coherent manybody ground state which 
may be closely related to the Abrikosov's vortex solution. Feynman's vortex solution of the 
superfluid He-4 may be a good starting point. More details will be published elsewhere. 26 
Gor'kov's formalism may also be generalized to obtain a coherent superposition of the dif- 
ferent types of the pairing. 17 ' 18 
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